arXiv:math/0512370vl [math.AG] 15 Dec 2005 


Elementary proof of the B. and M. Shapiro 
conjecture for rational functions 

Alex Eremenko and Andrei Gabrielov* 

October 10, 2005 


Abstract 

We give a new elementary proof of the following theorem: if all 
critical points of a rational function g belong to the real line then 
there exists a fractional linear transformation cj) such that cf) o g is 
a real rational function. Then we interpret the result in terms of 
Fuchsian differential equations whose general solution is a polynomial 
and in terms of electrostatics. 

One of the many equivalent formulations of the Shapiro conjecture is the 
following. Let f = (/i,..., /p) be a vector of polynomials in one complex 
variable, and assume that the Wronski determinant lE(f) = W{fi,..., fp) 
has only real roots. Then there exists a matrix A G GL(p, C) such that iA 
is a vector of real polynomials. 

This conjecture plays an important role in real enumerative geometry 
HSl 113, theory of real algebraic curves [S] and has applications to control 
theory cnnn. There is a substantial numerical evidence uni in favor of the 
conjecture. 

In pp we proved the Shapiro conjecture in the first non-trivial case p = 2. 
The proof was quite complicated, and its main drawback from the point of 
view of generalizations to higher dimensions was the use of the Uniformiza- 
tion theorem. 

In this paper we give a new proof, not using the Uniformization theorem. 
The new proof is also much simpler than the arguments in p. 

*Both authors were supported by NSF. 
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Consider a non-constant rational fnnction g = / 1 // 2 , and assnme that 
the polynomials /i and /2 are co-prime. Then the degree of g is given hj d = 
max{deg /i, deg / 2 }, and the roots of the Wronski determinant W (f) = / 1/2 — 
/ 1/2 coincide with the critical points of g. Let ns call two rational fnnctions 
gi and g 2 eqnivalent if gi = (j) o g 2 for some fractional-linear transformation 
(j). Evidently, eqnivalent rational fnnctions have the same critical points. So 
onr resnlt is 


Theorem 1. If all critical points of a rational function are real then it is 
equivalent to a real rational function. 


It is enongh to prove this theorem for rational fnnctions with simple 
critical points. The general case then follows by a limiting process. 

It is known jB] that for given 2d —2 points in the complex plane in general 
position, there exist 


1 / 2d - 2 \ 

“'"“dV J’ 


( 1 ) 


the d-th Catalan nnmber of classes of rational fnnctions of degree d with these 
critical points. It tnrns ont that the general position assnmption in this resnlt 
can be removed if the given points are real. Moreover, the following resnlt 
tnrns ont to be eqnivalent to Theorem 1: 


Theorem 2. For any given 2d —2 distinct points on the real line, there exist 
exactly Ud distinct classes of rational functions of degree d with these critical 
points. 


It follows from Theorem 1 that each of these Ud classes contains a real 
fnnction. The assnmption that the critical points are real is essential in 
Theorem 2: for 2d — 2 complex points, the nnmber of rational fnnctions of 
degree d with these critical points can be less than Ud- 

Eqnivalence of theorems 1 and 2 was known for some time, see, for ex¬ 
ample, ng. 

To state a generalization of Theorem 2 to the case of mnltiple critical 
points, we recall the dehnition of Kostka numbers. Let a = (oi,..., Og) be a 
vector of integers satisfying 


Q 

l<aj <d-l, ^ Oj = 2d - 2. (2) 

i=i 


Consider the Yonng diagrams of shape 2 x (d — 1). They consist of two rows 
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of length d—1. A semi-standard Young tableau SSYT of shape 2 x (c?— 1) is a 
hlling of such a diagram by positive integers, such that an integer k appears 
Qk times, the entries are strictly increasing in columns and non-decreasing in 
rows. The corresponding Kostka number Kg, is the number of such SSYT. 
The number Kg does not change if the coordinates of a are permuted na 
Thm. 7.10.2]. 

Theorem 3. For given a satisfying m, and given real points Xi < X 2 < 

... < Xq, there are exactly Kg classes of rational functions of degree d with 
critical points at Xj of multiplicity aj. 

We obtain Theorem 2 as a special case when q = 2d — 2 and all aj = 1. 
Theorem 3 is true for generic complex Xj] this was derived by Scherbak ffH 
from the results in m- Theorem 3 was hrst proved in where a result 
from pp was used. We include a proof here to show that it can be achieved 
with the same elementary tools as theorems 1 and 2, and no heavy machinery 
from pp is needed. 

A self-contained proof of Theorems 1 and 2 is given in Section 1. In Sec¬ 
tion 2, we discuss two equivalent reformulations of these theorems, which in 
our opinion are of independent interest, and also might be helpful for proving 
the Shapiro Conjecture in higher dimensions. Many other reformulations of 
the Shapiro conjecture are contained in usi uni uni- 

We thank Boris Shapiro, Frank Sottile and Alexander Varchenko for their 
useful comments on this paper. 


1 Proof of the main theorems 

1.1 Wronski map 

We recall the necessary facts and dehnitions. Let f = (/i,/ 2 ) be a pair 
of linearly independent polynomials of degree at most d. They span a 2- 
dimensional subspace in the space of all polynomials of degree at most d, 
and thus dehne a point in the Grassmannian G = G{2, d 1). Two pairs of 
polynomials are equivalent if they span the same subspace. 

The Wronski determinants of equivalent pairs are proportional non-zero 
polynomials of degree at most 2d — 2. Classes of proportionality of such 
polynomials form a space Poly^'^”^ which can be identihed with the projective 
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space So taking Wronski determinant defines a map 

W : G(2, d+1)^ Poly^'^-^ (3) 

which is called the Wronski map. The real Grassmannian Gr or the real 
projective space PolyR consist of those points whose coordinates (coefficients 
of the polynomials) can be chosen real. It is clear that W sends Gr to 
Poly^”^. The Wronski map is a hnite regular map of compact algebraic 
manifolds, and its degree can be dehned as the number of preimages of a 
generic point. This number turns out to be the Catalan number Ud, see, for 
example, jB]. If the Grassmannian G is embedded in a projective space by 
the Plucker embedding, then the Wronski map W becomes a restriction of a 
linear projection on G. Thus the degree of W is the same as the degree of the 
Grassmann variety, that is the number of intersections of a generic subspace 
of codimension 2d — 2 with the Plucker embedding of G. 

Using this notation, theorems 1 and 2 can be restated as follows: 

1. The full preimage of a polynomial with all real roots under the Wronski 
map consists of real points in G. 

2. Every polynomial in Poly^'^”^ with distinct real roots has exactly Ud distinct 
preimages under the Wronski map. 

If a pair (/i, f 2 ) represents a point of G, then S' = / 1//2 is a non-constant 
rational function of degree at most d. If two pairs of polynomials represent 
the same point of G, then the corresponding rational functions are equivalent. 
This dehnes a map r from G into the set Rat'^ of equivalence classes of non¬ 
constant rational functions of degree at most d. This map is not injective 
because polynomials in a pair can have a common factor. More precisely, let 
Zq G G he the locus of points corresponding to pairs of polynomials having a 
non-constant common factor, and Zi the locus corresponding to polynomials 
of degree less than d. Put Z = ZqU Zi. Then 

r : G\Zo Rat'^ (4) 

is a bijection. The standard topology on G can be dehned as induced by 
the Plucker embedding, and the topology on Rat^ is of uniform convergence 
with respect to the spherical metric. The map r is continuous on G\Z but 
not continuous on the whole G. The following weaker continuity property of 
r holds on the whole Grassmannian. 


4 


Proposition 1. Let {pj) be a converging sequence in G, pj G G\Z and 
p = \mipj is represented by a pair of polynomials with a common factor q. 
Let zi,...,Zk be the roots of q. Then r{pj) r{p) uniformly on compact 
subsets of C\{2:i, .. ., Zk, cxd}. One does not have to include oo if the degree 
of p is the same as that of pj. 

The elementary proof is left to the reader. 

1.2 Nets of rational functions 

Let be the class of real non-constant rational functions g of degree at most 
d whose all critical points are real. Consider the full preimage 7 = 
where R = RU {cx)}. This preimage consists of simple analytic arcs which 
meet only at critical points. These arcs dehne a cell decomposition G{g) 
of the Riemann sphere whose 2-cells (faces) are components of P^\ 7 , 1- 
cells (edges) are components of 7 \{critical points} and 0 -cells (vertices) are 
the critical points. We choose some vertex uq and call it the distinguished 
vertex of G{g). The union 7 of edges and vertices is the 1 -skeleton of the 
cell decomposition. Such a cell decomposition G = G{g) has the following 
properties: 

(i) the 1 -skeleton of G contains R. 

(ii) G is symmetric with respect to R, 

(hi) all vertices belong to R and an even number of edges meet at each vertex. 

If C* = G{g) the even number in (hi) is twice the local degree of g at 
the critical point. Another important property of our cell decomposition is 
that the closure of every cell is homeomorphic to a closed ball of the same 
dimension. It follows that 

(iv) no edge can begin and end at the same point. 

Two cell decompositions Gi and G 2 satisfying (i)-(iv) and having distin¬ 
guished vertices Vq and Vq will be called equivalent if there exists a homeomor- 
phism 0 : P^ —P^ commuting with reflection with respect to R, preserving 
orientations of R and P^, mapping cells of Gi onto cells of G 2 and Vq to Vq. 
An equivalence class of cell decompositions satisfying (i)-(iv) will be called 
a net. The number of faces of a net is even, we denote it by 2d and call the 
positive integer d the degree of the net. If C = G(g) then degC = deggf. A 


5 



net of degree d has 2d — 2 edges disjoint from the real axis. Using the Uni- 
formization theorem, we proved in DQ that each net comes from a rational 
function, and the critical points of this rational function can be arbitrarily 
prescribed, but we do not use this result here, and in fact it will be deduced 
from our theorems 1 and 2 in the end of Section 1. 

We need the following elementary 

Proposition 2. Let {pj) be a convergent sequence in Gr,, and p = limp^. 
Let Qj = r{pj) be the corresponding sequence of rational functions. Then the 
sets 7 j = g~^{lV) converge in the Hausdorff metric to the set 7 = 5 f“^(R), 
where g = r{p). 

This is a simple corollary of Proposition 1, and the details are left to the 
reader. 

Corollary 1 . Suppose that gt : t ^ [0,1] is a continuous path in and each 
gt has 2d — 2 simple critical points. Let one of these critical points be Vo{t), 
a continuous function oft. Then the net of gt with distinguished point Vo{t) 
is independent oft. 

Corollary 2. Let pt be a continuous path in the Grassmannian G, paramet¬ 
rized by [0,1]. Suppose that gt = r{pt) belong to and have critical points 
Xo(t),... ,Xn(t), such that Xj(t) 7 ^ Xi(t) for 0 < j < i < n and 0 < t < 1, 
while for t = 1 we have a:o(l) = ®i(l) ^j(f) 7 ^ ^i(f) fof 1 < J < * < 

Then the degrees of gt are equal for t G [0,1), and the degree of gi is less 
than the degree of gt, t G [0,1) if and only if the net of go contains an edge 
from Xq to Xi. 

Proof. According to Proposition 2, the sets jt = vary continuously 

in the Hausdorff metric. If there is an edge connecting Xo{t) and Xi{t), the 
limit of this edge as t > 1 cannot be a loop because of the property (iv) of 
the nets. So the limit belongs to R and thus the limit cell decomposition has 
fewer faces than the cell decomposition C{gi). 

In the opposite direction, if this limit has fewer faces, some edge has to 
disappear in the limit, and this can only be an edge from xq to Xi. □ 

1.3 Thorns 

Let C Poly^“^ be the subset consisting of polynomials whose all roots 

are real. Then has non-empty interior. Here we construct an open 
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subset of such that for every polynomial p in this subset, the full 

preimage W~^{j)) consists of Ud distinct real points in G. The existence of 
such a subset was established by Sottile [Hj, but we give a more precise 
description of this set following 

We £x d > 2. Consider pairs of integers 


0 < fci < ^2 < d 


(5) 


and pairs q = (gi, ^ 2 ) of real polynomials 

qi{z) = + ... + 

q2{z) = 0‘2,d-G'^ ^ + • • • + 


( 6 ) 


Suppose that all coefficients aij are strictly positive, all roots of the Wron- 
skian determinant W{q) = W{qi,q2) belong to the semi-open interval 
(—1,0] C R, and those roots on the open interval (—1,0) are simple. The 
set of all such polynomial pairs m will be denoted by b{ki, ^ 2 ). The greatest 
common factor of {^ 1 ,^ 2 } is z’^^. 

It is easy to see that the representation of a point of Gr by a pair in 
b{ki,k 2 ) is unique. Setting ki = 0 and A ;2 = 1 we obtain an open subset 
6(0,1) of Gr. 

The multiplicity of the root of W (q) at 0 is /c = /ci -|- ^2 — 1- We enumerate 
the negative roots of W (q) as 


- X2d-2 < -X2d-3 < ■ ■ ■ < -Xk+i < 0. (7) 

Let e be a positive increasing function on [0,1] satisfying e(0) = 0, e(x) < 
X, for X G (0,1]. The set of all such functions will be denoted by E. A thorn 
T of dimension n is a region in of the form 

T{n, e) = {{yi, 0 <?/„ < e(l), 0 < yk < e{yk+i), 1 <k <n-l}. 

Let w{k,T) be the set of polynomials of the form 

p{x) = x^{x + X2d-2){x -t- X2d-3) . . . (x Xfc+l), ( 8 ) 

where the vectors {x 2 d- 2 , ■ ■ ■, Xk+i) belong to a thorn T of dimension 2d—2—k. 
The set 6(d — 1, d) consists of a single pair qi = z'^~^, q 2 = z^. 

Consider the following two operations F\ i = 1, 2. For each pair (gi, g 2 ), 
of the form © operation F* adds to the polynomial g^ one term az’^* 
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where a > 0 is a small parameter, and leaves the other polynomial of the 
pair nnchanged. So each operation F* increases the total nnmber of non-zero 
coefficients of a polynomial pair by one: 

Fl:b{k)^b{k-e.), f;(q)=q + o2*<-‘e„ ( 9 ) 

where k = (^ 1 ,^ 2 ) and ( 61 , 62 ) is the standard basis in R^, and a > 0 is a 
small parameter whose range may depend on q. The following rnle has to 
be observed: 

Rule. Operation F* is permitted on b{ki, ^ 2 ) if o,nd only if the outcome of 
this operation does not violate inequalities ©• 

In other words, operation F^ is permitted on h{ki,k 2 ) if ki > 0, and 
operation F^ is permitted on b{ki, ^ 2 ) if k2 > ki + 1 . 

The following result, which is a part of 0 Proposition 8 ] shows, among 
other things, that the F* are well dehned, that is their result indeed belongs 
to some b{kl, k^) for sufficiently small values of the parameter a. 

Proposition 3. Suppose that k = {ki,k 2 ) and i G {1,2} satisfy the Rule 
above. Suppose that for some thorn T of dimension 2d — 2 — k a set U C 
b{ki, /C 2 ) is given, such that the map q 1 —> IT(q) : —> w{k,T) is surjective. 

Then there exists a thorn T* of dimension 2d —1 — k and a setU* C 6(k — e*), 
such that every q* G U* has the form F*(q), q G 17, where F^ is defined in 
©, and a > 0, and the map 

q* :U* ^w{k-l,T*) (10) 

is surjective. 

Proof. We follow |Sl Section 2], First we state three elementary lemmas 
about thorns. 

Lemma 1. Intersection of any finite set of thorns of same dimension con¬ 
tains a thorn of the same dimension. 

Proof. Take the minimum of their dehning functions. □ 

Lemma 2. Let T = T{n, e) be a thorn of dimension n in R" = {(xi,..., Xn)}, 
and U its neighborhood in R”'+^ = {(xq, ... , x„)}. Then = 17fl R>o^ con¬ 
tains a thorn T{n -\- 1, ei). 

Proof. There exists a continuous function Sq : T ^ ^> 0 ; such that 17+ 
contains the set {(xo,x) : x G T, 0 < Xq < (5o(x)}. Let 6{t) be the minimum 


of 5o on fhe compact subset {x G T(n, e/2) Xi > t} of T. Then there 
exists eo E E with the property Cq < S. If we dehne ei = min{e/2,eo}, then 
T(n+ l,ei) C □ 

Lemma 3. Let T = T(n + 1, e) be a thorn of dimension n + 1, and 
h : T ^ ^-> 0 ^) (a^o,x) i—> {yo{xo,x),y{xo,x)), a continuous map with the 
properties: for every x such that (a;o,x) G T for some Xq > 0, the function 
Xq h->• yo{xQ, x) is increasing, and lima-g^o y( 2 ^ 0 ; x) = x. Then the image h(T) 
contains a thorn. 

Proof. We consider the region D G consisting of T, its reflection T' 
in the hyperplane xq = 0 and the interior with respect to this hyperplane of 
the common boundary of T and T'. The map h extends to T' by symmetry: 
h{—Xo,x) = —h{xo,x), (xo,x) G T, and then to the whole D by continuity. 
It is easy to see that the image of the extended map contains a neighborhood 
U of the intersection of D with the hyperplane xq = 0. This intersection is 
a thorn Ti in = {(a;o,x) G : xq = 0}. Applying Lemma 2 to this 

thorn Ti, we conclude that U~^ contains a thorn. □ 

We continue the proof of Proposition 3. 

Let us £x q G U, and put W = ITq. As hP G w{k, e), we have ord W = k, 
where ord denotes the multiplicity of a root at 0. Let cz^ be the term of the 
smallest degree in W{z). Then c > 0, because all roots of W are non-positive. 
In fact, 

c={k2- ki)a2^k2aiM > 0- (H) 

We fix i G {1,2} satisfying the Rule above and define W* = hPq*, where 
q* = i/)(q). Then ordhP* = k — 1 and the term of the smallest degree in 
W*{z) is c*z^~^, where 

c* = a{kl - A;{)a3-i,fc3_, > 0, (12) 

We conclude that when a is small enough (depending on q), the Wronskian 
W* has one simple root in a neighborhood of each negative root of W, and 
in addition, one simple negative root close to zero, and a root of multiplicity 
/c — 1 at 0. To make this more precise, we denote the negative roots of W 
and hP* by 

- X 2 d -2 < ... < -Xk+i and - y 2 d -2 < ... < -yk+i < -Vk, (13) 
where yj = yj(a). We have 

yj{0) = Xj, for I < j <n, and 1/^(0) = 0. (14) 
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Furthermore, if a is small enough (depending on q) 

a H->• yk{a) is increasing and continuous. (15) 

The set w{k,e) is parametrized by a thorn T = T{2d — 2 — k,e), where 
X = (xfc+i,..., X 2 d- 2 )- There exists a continuous function 6 : T —!■ R>o, such 
that 

q* e 5(k*), for aG(0, (5(x)), x G T. (16) 

It remains to achieve m by modifying the thorn T. Consider the set 

U* = {q* = F„(q,) : X G T, a e (0, 5(x))} C 6 (k*), (17) 

where qx G 17 is some preimage under W of the polynomial (jH)) with 


{xk+i, ■ ■ ■, X 2 d- 2 ) = X. Such preimage exists by assumption of Proposition 3 
that the map q h-> IFq, U w{k, T) is surjective. We apply Lemma 2 to the 
half-neighborhood (ED of T, with Xk = a, to obtain a thorn Ti{2d — k — l, ei). 
Then we apply Lemma 3 to the map h : Ti ^ dehned by yj = 

yj{xo,x.), where yj are as in (fT!?|l . and Xk = a. 

This map h satishes all conditions of Lemma 3 in view of ED and (USD. 
This proves (ED and Proposition 3. □ 

We begin with the single element of h{d — 1, d) and apply operations F* 
in some sequence, obeying the Rule above, while possible. As every step 
decreases A; by 1, the total number of steps will be 2d — 2. We describe the 
sequence of steps by a sequence cr of I’s and 2’s of length 2d —2. The number 
i on the n-th place in this sequence indicates that operator F* was applied 
on n-th step. The Rule above translates to the following characterization of 
all possible sequences a: 

a) The numbers of I’s and 2’s in a are equal. 

b) In each initial segment of a the number of I’s is not less than the 
number of 2’s. 

Such sequences are called ballot sequences (for two candidates), see, for 
example, in|. The number of ballot sequences of length 2d — 2 is the Catalan 
number Ud- 

Applying 2d —2 times Proposition 3 according to each ballot sequence a 
we obtain in the end an open subset Ua- C 5(0,1) which is mapped surjectively 
by the Wronski map onto w{0,T) for some thorn T. As the intersection of 
any hnite set of thorns of the same dimension contains a thorn of the same 
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dimension by Lemma 1, we may assume that this set w(0, T) is the same for 
all sequences a. 

When applying Proposition 3 we can make the range of parameter a as 
small as desired; using this we can assure that the sequence of coefficients 
of the pair © is monotone: the coefficients decrease in the order of their 
appearance. This implies that the sets U„ with different a are disjoint. 

So we obtain Ud disjoint open sets and each of them is mapped onto 
r(;(0, T) continuously and surjectively by the Wronski map W. As the number 
of preimages of any point under W is at most Ud we conclude that the maps 
W : Ua ^ w{0,T) are homeomorphisms for all a. 

Thus each point of the open subset tn(0, T) C has exactly Ud preim¬ 

ages under W and all these preimages are real. Each of these preimages cor¬ 
responds to an analytic branch of the inverse W~^ on tc(0, T). The branches 
are enumerated by ballot sequences. 

1.4 Completion of the proof 

Let us £x a thorn T such that each polynomial in w{0,T) has Ud different 
real preimages under the Wronski map, as in the end of the previous section. 

To each of these preimages q = {qi,q 2 ) corresponds a rational function 
r(q) = qi/q 2 in with 2d — 2 distinct real critical points which has a net 
7 (q). We take the rightmost critical point of these functions as distinguished 
vertices of the nets. 

We claim that all these Ud nets are different. To prove the claim we just 
show how to determine the net from the ballot sequence and vice versa. 

Proposition 4. Let k = 0, and let p be a polynomial in w(0, T) of the form 
(0- Let q = {qi,q 2 ) be a polynomial pair as in 0 corresponding to a point 
in W~^{p), with the ballot sequence a, and g = r{q) = qi/q 2 . 

Then the net of g contains an edge between Xm and Xm+i if and only if 
the m-th member of the sequence a is 1. 

Proof. It is enough to investigate what happens to a net when an operator 
pt Qf Proposition 3 is applied. We see from © and o that the degree of 
qi/q 2 increases if and only if i = 1. Corollary 2 of Proposition 2 says that 
this happens if and only if the net has an edge between Xm and Xm+i- This 
proves Proposition 4. □ 

So we obtained a polynomial po £ of degree 2d — 2 with 2d — 2 
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real roots whose preimage under the Wronski map consists of Ud pairs with 
different nets. Each of these preimages corresponds to a holomorphic inverse 
branch of the Wronski map in w{0,T). Let pi be any real polynomial of 
degree 2d — 2 with 2d — 2 real roots. Then there exists a path pt : t & [0,1] 
in connecting po and pi, such that all pt are polynomials with 2d — 2 

distinct roots. For example one can connect the corresponding roots of po 
and Pi linearly. We do analytic continuation of all the inverse branches of 
the Wronski map along this path. As critical points of our rational functions 
cannot collide (because the zeros of their Wronskians pt do not collide), their 
nets do not change under the continuation. Suppose that this analytic con¬ 
tinuation to t = 1 is impossible. Let to be the smallest singular point. Then 
PtQ is a ramified value of the Wronski map and the full preimage W~^{ptQ) 
consists of fewer than Ud points. This full preimage still consists of real 
rational functions with all critical points real and distinct, so the nets are 
dehned for all elements of this preimage. This means that at least two one- 
parametric families of rational functions with different nets tend to the same 
function with 2d —2 distinct critical points, which is impossible by Corollary 
1 of Proposition 2. 

This proves theorems 1 and 2. 

This proof clearly implies that for any given net there exists a unique class 
of real rational functions with all critical points real, and the critical points 
of these functions can be chosen arbitrarily, the result which was established 
in pp with the help of the Uniformization theorem and rather complicated 
topological considerations. 

To prove Theorem 3, we notice that 1-skeleton of every net of degree d 
can be obtained as the limit of 1-skeletons of nets with 2d — 2 vertices. So 
every net of degree d actually occurs as a net of a real rational function of 
the class of degree d with all critical points real. Counting the nets of degree 
d with q vertices of degrees 2ai..., 2aq gives the Kostka number (see, for 
example, 0 Lemma 3]. So there are at least classes of rational functions 
of degree d with prescribed real critical points. On the other hand, Schubert 
calculus CH shows that there are at most classes of rational functions 
with any prescribed critical points of multiplicities Oi,... , 0 ^. This proves 
Theorem 3. 

Corollary. To each net of degree d corresponds exactly one class of real 
rational functions of degree d with prescribed real critical points. 
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2 Fuchsian equations 


2.1 Equations with all polynomial solutions 

Suppose that distinct points ai,... ,an in the complex plane are given. We 
want to describe the set of equivalence classes of polynomial pairs {yi,y 2 ), 
such that 

W{yi,y2) {z - ai).. .{z - a„). 

This time we do not specify in advance the degree of polynomials yi,y 2 but 
it is easy to see that it is at most n + 1. As fT(|/i, 1 / 2 ) has only simple zeros, 
the polynomials yi and y 2 are co-prime. 

To approach this problem, we introduce a special parametrization of 
equivalence classes of polynomial pairs whose Wronskian has prescribed ze¬ 
ros. Recall that for any two linearly independent fnnctions yi and y 2 one can 
write a second order linear differential eqnation which has these two fnnctions 
as solntions: 

y Vi y2 
y' y'l 1/2 
y” y'l y'i 

where 

A = W(|/i,|/ 2 ), B = -A' and C = W{y[,y'.,). 

If yi and y 2 are polynomials, then A, B, C are also polynomials, and we 
have 

degR < deg A — 1, degC<degA —2. (19) 

These conditions are eqnivalent to regnlarity of the singular point at inhnity. 
Introdncing two rational fnnctions P = B/A and Q = C/A, we conclnde 
that P{oo) = Q{oo) = 0. Furthermore, P = —A'/A, and A has only simple 
zeros, so all residnes of P in C are eqnal to —1. Denoting the residnes of Q 
by Xj, we obtain 


= + Bj/'+ Cj/= 0, (18) 


P{z) 


E 

i=i 


and Q{z) = 


Xn 


j=i 


( 20 ) 


Evidently, P, Q and Xj depend only on the eqnivalence class of the polynomial 
pair {yi,y 2 ). 
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Now we write the conditions on Xj which express the fact that all solutions 
of the differential equation 


y" + Py' + Qy = o 


( 21 ) 


are polynomials. 

In a neighborhood of a singular point a^, our equation can be written in 
the form: 

{z - ak)y'' + Pk{,z)y' + Qk{,z)y = 0, (22) 

with 

Pk{z) = -1 + Pk{z - ttk) + 0{z - ttkf, where —-—, 

and 

Qk{z) = Xk + qk{z - ttk) +0{z - ttk)^, where ^—. 

T^k 

Our differential equation (1221) has two linearly independent polynomial solu¬ 
tions without a common factor, so it has a polynomial solution of the form 

y{z) = 1 -h ci{z - ak) + C2{z - akf + 0{z - akY- 

Differentiating this, we obtain 

y\z) = Cl -1- 2c2{z - ttk) + 0 {z - ttkY, 

y"{z) = 2 c2 + 0 {z- ak). 

Substituting this to (j221), we obtain 


Cl Xj^ , 


and 


PkCi + XkCi + qk = 0, 


so 

= -PkXk - qk- 


(23) 
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Recalling the expressions for PkiQki we obtain the following necessary con¬ 
dition for the equation ( 1221 ) to have two linearly independent polynomial 
solutions: 

xl = ^——= (24) 

Proposition 5. Condition is necessary and sufficient for a differential 
equation dH with coefficients (US) to have two linearly independent polyno¬ 
mial solutions. 

Proof It remains to prove sufficiency. From (I22|l we conclude that all 
singular points in C are regular, with exponents 0 and 2. Condition (El 
guarantees that there is a power series solution corresponding to the smaller 
exponent. This implies that there are two linearly independent holomor- 
phic solutions in a neighborhood of each singular point. Thus all solutions 
are entire functions. By a theorem of Halphen which can be found in [71 
15.5], if the general solution of the equation ([H, where A,B,C satisfy 
degR < deg A, degC < deg A, is a meromorphic function in C, then this 
meromorphic function has to be of the form 


where Rj are rational functions. But asymptotics at inhnity shows that under 
the stronger condition dH on the degrees of the coefficients, the exponentials 
cannot be present. So the general solution is an entire rational function, that 
is a polynomial. □ 


Proposition 6. Every solution of the system {"EJ) has the following proper¬ 
ties: 


J2xk = 0, (25) 

k=l 


{n +if - = s'^, (26) 

k=l 

where s is an integer such that n -\- s is odd, and 1 < s < n -|- 1. 


This integer s is the local degree at inhnity of the rational function which 
is the ratio of two linearly independent solutions of EH) with P and Q as in 

(EOl)- 


15 




Proof of Proposition 6. By Proposition 5, for every solution xi,..., of 
(El, all solutions of the differential equation El are polynomials. Suppose 
that 


P{z)=pz ^ + 0{z ^), Q{z) = qz ^ + q*z + 0{z z ^ oo. 

Then 

n n 

p=-n, q = '^Xk and q* = y^^Xkttk- 

k=l k=l 

Substituting into El a polynomial 

y{z) = z^ + ..., 

we hrst obtain 

qz'^~^ + 0{z‘^~‘^) = 0, z —> cx), 

so g = 0 which proves (El- Then we obtain 

d{d - l)z‘^-^ + dpz‘^-^ + q*z^-^ + 0{z^-^) = 0, 
so 

— d{n + 1) + g* = 0. 

This equation has two solutions, the possible degrees of polynomials: 

n + 1 ± A/(n + 1)^ — 4g* 
dl,2 = -^-. 

As all solutions of our differential equation are polynomials, there are two 
solutions of different degrees. This implies that both di and ^2 are integers 
which implies fj2t)j) . Notice that di — d 2 = s. 

Rational function yi/y 2 is locally s to 1 at inhnity, from which the in¬ 
equality 1 < s < n -|- 1 follows. □ 


So our theorems 1 and 2 are can be restated as 

Theorem 4. If all Ok are real then all solutions of are real. More¬ 
over, each solution (xi,... ,x„) can be analytically continued as a function of 
parameters (oi,..., On) in the region oi < 02 < ... < a„. 

Solutions (xi,...,x„) correspond to classes of rational functions of all 
degrees d G [n/2 -|- 1,?7, -|- 1] if n is even and d G [(n -|- l)/2 -f 1,?7, -|- 1] if 
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n is odd, having simple critical points at Oi,... ,a„ and possibly a critical 
point at infinity. The order of the critical point at infinity is s — 1 where s is 
defined in 

System (1^ has a trivial solution xi = ... = = 0 which corresponds to 

the polynomial of degree n + 1 with critical points oi,..., a„. The opposite 
case is that s = 1 (then n is even) and we have rational functions with n 
prescribed simple critical points. This case is characterized by the condition 


q* = (n^ + 2n)/4. 


System can be easily generalized to p-tuples of linearly independent 
polynomials satisfying a Fuchsian differential equation of order p. An analog 
of Theorem 4 for this case would be equivalent to the Shapiro conjecture for 
p polynomials. For example, for p = 3 one obtains the following system of 
equations with respect to Xk and Uk- 




2.2 Equilibria of electric charges in the plane 

Following Stieltjes we state an extremal problem of potential theory 
which is equivalent to (ED). 

Suppose as above that 


A(z) = {z — ai) ■■■ {z — a„), and B = —A' 


(27) 


in m, and let us look for polynomials C such that m has a polynomial 
solution 


y{z) = {z - Zi) .. .{z - Zm) with Zj G C\{ai,..., Un}. (28) 
We substitute in (fTK|l and put z = Zk- 

A{zk)y"{zk) - A'{zk)y{zk) = 0 , 


or 


y''{zk)/y'{zk) - A'{zk)/A{zk) = 0 , 
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which is equivalent to 




1 


Zk - Zj 



1 




0 . 


(29) 


Equation (Ei has the following physical interpretation. Positive unit charges 
+1 are hxed at the points a^, and negatively charged particles of charge 
—2 each at the points Zk are allowed to move in the plane. The particles 
are repelled or attracted according to their charges, and the force is inverse 
proportional to the distance^. So the force between z and w is c/{z — w), 
where the bar stands for the complex conjugation and c is a real constant 
depending on the charges. Then the equilibrium condition is expressed by 
the system of equations (P|) . The energy of such conhguration is 

^ ^ iQg Uij,ky.j^k\^k - zj\^ ^ 

n(i,fc) \^k — Oji 

The function U{z,sl) under the logarithm is called the master function in 
[13 El im US- Here z = {zi,, Zm) is the variable and a = (oi,..., a„) the 
parameter. Equilibrium conhgurations are the critical points of the master 
function. It is evident from physical considerations and easy to prove that 
all these equilibria are unstable. 

So for a given polynomial A as in (EH), an equation of the form 


A-i/' - A'-i/ + Cy = 0 


(30) 


has a non-trivial polynomial solution y if and only if 

C={-Ay" + Ay')/y, 

with y given by (OHll and {zi,..., Zm) satisfying (ED). 

Lemma 4. If the equation PI has one non-trivial polynomial solution, then 
all its solutions are polynomials. 

This is [TTIl Lemma 7]. We include a proof for the reader’s convenience. 

'^The particles in the plane can be imagined as infinite uniformly charged wires per¬ 
pendicular to the complex plane and interacting by the Coulomb law. 
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Proof. We apply the usual method of hnding a second linearly indepen¬ 
dent solution when one, say yi, is given. We obtain 



(31) 


This can have ramihcation points only at the zeros of yi. But as a solution of 
a differential equation dSOD, r /2 can only have ramihcation points at the zeros 
of A. As these two sets are disjoint, y 2 is an entire function. Now, using 
the expression m we conclude that y 2 is a polynomial. So all solutions are 


polynomials. 


□ 


Remark. An alternative way to derive (I29|l is to write the condition that 
all residues of the integral in m are equal to zero. 

We conclude that some solutions of (Ei are not isolated; they occur in 
complex 1-dimensional families corresponding to conhgurations of zeros of 
all polynomial solutions of dsni). However, each two-dimensional space of 
polynomials contains a one-dimensional subspace consisting of polynomials 
of smaller degree than generic polynomials in this space. Such polynomial 
of the smallest degree gives an isolated solution of (f^ . If n = 2d — 2 
and m < d — 1, we see from m that deg 1/2 > deg|/i, and we conclude 
that the critical points of the master function E are all isolated in this case. 
These facts about critical points of the master function and their relation to 
Fuchsian equations with polynomial solutions were discovered in They 
can be generalized to Fuchsian equation of arbitrary order. The equation 
(Ei is a special case of the “Bethe Ansatz equation” in na. 

Suppose now that EDI) has a polynomial solution. Then all solutions 
are polynomials and the ratio of any two linearly independent solutions is 
a rational function with all critical points real. Such a rational function is 
equivalent to a real rational function by Theorem 1, so we can hnd two real 
linearly independent solutions. Then C* is a real polynomial, and the solution 
of smallest degree is proportional to a real polynomial. Thus Theorem 1 has 
the following consequence: 

Theorem 5. For given real points Oi,..., a„, each isolated equilibrium con¬ 
figuration zi,... ,Zm in (113 is symmetric with respect to the real line. 

It is also easy to deduce Theorem 1 from Theorem 5. 
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